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We summarise recent results on the properties of gluons, quarks and light mesons from the Green's functions 
approach to QCD. We discuss a self-consistent, infrared power law solution for the Schwinger-Dyson equations of 
the IPI-Greens functions of Yang-Mills theory. The corresponding running coupling has a universal fixed point at 
zero momentum. Based on these analytical results a truncation scheme for the coupled system of Schwinger-Dyson 
equations for the propagators of QCD and the Bethe-Salpeter equation for light mesons has been formulated. We 
compare numerical results for charge eigenstate vector and pseudoscalar meson observables with corresponding 
lattice data. The effects of unquenching the system are found to be small but not negligible. 



One of the most fascinating problems of QGD 
is to find a low energy description of colourless 
bound states (hadrons) in terms of their non- 
perturbative, coloured constituents (quarks and 
gluons). Lattice simulations are not entirely sat- 
isfactory in this respect. They provide values 
for the global properties of hadrons (masses, de- 
cay widths etc.,), but they may not be capable 
to provide enough information on their internal 
structure to understand all dynamical aspects of 
low energy QGD. An alternative field theoretical 
and relativistic method which is well suited to de- 
liver this information is the Schwinger-Dyson and 
Bethe-Salpeter formalism |1I2| . In principle this 
approach allows one to derive meson properties 
directly from the fundamental building blocks of 
the field theory, the Green's functions. Lattice 
simulations and the Green's functions approach 
are complementary to each other in several re- 
spects. Lattice simulations are ab initio whereas 
the Green's functions approach has to rely on an 
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(educated) approximation scheme. The Green's 
functions approach is continuum based. It allows 
for analytical investigations in the infrared and 
all aspects of chiral symmetry and its breaking 
pattern are respected. In this talk we summarise 
recent results in the Green's functions framework 
that dwell on these advantages. We discuss an- 
alytical results on the infrared exponents of the 
IPLGreen's functions of SU(Afc)-Yang-Mills the- 
ory in Landau gauge. We report on numerical re- 
sults for the ghost, gluon and quark propagators 
as well as light meson observables in a truncation 
scheme that is based on the analytical findings. 
Our focus in particular is on unquenching effects 
due to light quark loops in the gluon polarisation. 

1. Infrared exponents for the Green's func- 
tions of Yang-Mills theory 

The infrared behaviour of the Green's functions 
of SU(Arc)-Yang-Mills theory is related to confine- 
ment in several ways. A particularly interesting 
example is the Kugo-Ojima confinement criterion 
of a well-defined global colour charge. The crite- 
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rion is satisfied in Landau gauge if tfie dressing 
function of tlie ghost propagator is singular and 
the gluon propagator is finite or vanishes in the 
infrared. Provided BRST-symmetry is conserved 
nonperturbatively, the cohomology of the BRST- 
operator then defines a physical state space with 
colourless states only (51. 

A convenient starting point to investigate 
the infrared behaviour of one-particle-irreducible 
(IPI) Green's functions is the Schwinger-Dyson 
equation (SDE) for the ghost-gluon vertex, shown 
diagrammatically fig.G] The dressed ghost-gluon 
vertex r°''° = r^(p, q)f°''"' can be represented by 
the two tensor structures 

T^ip,q) ^ p^Aip-'^q^) + k^B{p'',q^), (1) 

where p^ and are the momenta of the out- 
going and incoming ghost and fc^ = — — 
is the gluon momentum. In Landau gauge, the 
momentum q^ of the incoming ghost factoriscs 
from the vertex dressing, as can be seen from 
fig. ^ Since the gluon propagator 1?^,^ is trans- 
verse in Landau gauge, its contraction with the 
bare ghost-gluon vertex in the loop of the SDE 
gives lfj,Df^^{l - q) ^ q^LD^,y{l - q). 
Let us assume for the moment that 
(I) the loop-integral is finite in the infrared. 
(We come back to this assumption in the para- 
graph below eq. ((HJ.) We then observe that the 
dressing of the full ghost-gluon vertex vanishes if 
q^ goes to zero, 

T,,{p,q)^p^{A-B)-q^B p^, (2) 

and thus neither A{p'^,q'^) nor B{p^,q^) can be 
singular in this limit Since factorisation of 
the other ghost momentum is trivial, the same 
happens for — > 0. One thus concludes that 
the dressing of the ghost-gluon vertex is finite in 
the infrared and may be well approximated at 
small momenta by the bare vertex. This has been 
confirmed by lattice and SDE-studies )5I6I7| . 

A finite ghost-gluon vertex at small momenta 
admits the following power law solution for the 
ghost-SDE in the infrared: Writing the ghost and 
gluon propagators as 




Figure 1. Schwinger-Dyson equation for the 
ghost-gluon vertex. 




Figure 2. Schwinger-Dyson equation for the ghost 
propagator. 

. (,„„-£*)5i|!). (3) 

one finds power laws for the ghost and gluon 
dressing functions with interrelated exponents 
given by 

[NB: This can be checked easily by just counting 
anomalous dimensions on both sides of the equa- 
tion. Note that the loop-integral is dominated by 
momenta of the same magnitude as the external 
momentum. Thus, for small external momenta 
one can replace the propagators in the loop by 
their infrared approximation, eq. .] In this no- 
tation the Kugo-Ojima criterion translates to the 
condition k > for the ghost dressing function 
and K > 0.5 for the gluon dressing function. On 
general grounds, the exponent k is known to be 
positive 0] , independent of any truncation of the 
SDEs. For a bare ghost-gluon vertex in the in- 
frared one obtains k = (93 - yT20T)/98 « 0.595 
|8l9j . which satisfies both criteria. This specific 
value of K is found to vary only slightly for a large 
class of possible dressings of the ghost-gluon- 
vertex 0. Similar values have been determined 
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Figure 3. Exact Schwinger-Dyson equation for 
the three-gluon vertex and lowest order in a skele- 
ton expansion of the four- and five-point func- 
tions. All internal propagators in the diagrams 
are to be understood as fully dressed. 



from exact renormalisation group equations |1U! . 

The basic idea to determine the infrared be- 
haviour of the other (IPI) Green's functions is to 
investigate their Schwinger-Dyson equations or- 
der by order in a skeleton expansion (i.e. a loop 
expansion using full propagators and vertices). 
This program has been carried out in ref. JT] . It 
turns out that in this expansion the Green's func- 
tions can only be infrared singular, if all external 
scales go to zero. Thus to determine the degree 
of possible singularities it is sufficient to investi- 
gate the SDEs in the presence of only one external 
scale <C ^QCD- example we summarise 

the treatment of the SDE for the three-gluon ver- 
tex. In fig. O we show the full equation as well 
as the lowest order in a skeleton expansion of the 
four-and five-point functions. In the presence of 
one (small) external scale the approximated SDE 
has a selfconsistent power law solution given by 

r'^ip') ~ (pT'^- (5) 



[NB: Again this can be seen easily by counting 
anomalous dimensions on both sides of the equa- 
tions. The leading diagram on the right hand side 
is the one involving ghosts, diagram (a), the oth- 
ers are less singular (recall k > 0). The loops are 
again dominated by momenta of the same magni- 
tude as the external scale.] One can see by induc- 
tion that this solution is also present if terms to 
arbitrary high order in the skeleton expansion are 
taken into account. Thus the skeleton expansion 
is stable wrt. the infrared solution of the SDEs. 
This technique can also be applied to any other 
SDE. A self-consistent solution of the whole tower 
of SDEs is then given by 

Here r"'™(p^) denotes the infrared leading dress- 
ing function of the IPI-Green's function with 2n 
external ghost legs and m external gluon legs. By 
counting anomalous dimensions it can be checked 
easily that the expression eq. © indeed solves the 
full three-gluon vertex SDE in fig. O selfconsis- 
tently. Furthermore, inserting r^'^{p^) ^ (p^)^" 
together with the power laws Q into the SDE for 
the ghost-gluon vertex, fig. ^ one can verify the 
assumption (I) that the loop-integral of the ver- 
tex dressing is indeed finite in the infrared. Thus 
cq. (O is a truly selfconsistent infrared solution of 
the tower of SDEs. [NB: It is worth mentioning 
that the solution jnj also has the correct scal- 
ing behaviour such that the Slavnov- Taylor iden- 
tities of the renormalisation constants are sat- 
isfied. Since the theory is multiplicative renor- 
malisable these functions scale with the renor- 
malisation point in the same way as the IPI- 
functions with the^ external scale p^. E.g. the 
relation Z1/Z3 — Zi/Z^ between the three-gluon 
vertex, gluon propagator, ghost-gluon vertex and 
ghost propagator renormalisation constant leads 
to Zi{^'^) = {^'^)~^'^, which agrees with eq. ©.] 
Certainly, selfconsistency is not enough to es- 
tablish eq. as the 'true' solution of Yang-Mills 
theory in the infrared, since there may be other 
selfconsistent solutions of the SDEs. However, 
one may argue that the solution ® has an in- 
teresting property that qualifies it as a promising 
candidate: it leads to qualitative universality of 
the running coupling in the infrared. Renormal- 
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isation group invariant couplings can be defined 
from either of the primitively divergent vertices 
of Yang-Mills-theory, i.e. from the ghost-gluon 
vertex (gh), the three-gluon vertex (3g) or the 
four-gluon vertex (Ag) via 

a'^p') ^ '^^r-^(j?)fZ\vY-"^^, (8) 
a'^ij?) = 1^ r-\v')f Z\f ) ^ . (9) 

Using the SDE-solution © it is easy to see that 
all three couplings approach a fixed point in 
the infrared. The constants Ci may be different 
for each coupling and depends on the respective 
choice of the tensor component used to extract 
the vertex dressing functions F (this ambiguity is 
well know in the literature ) ■ For the coupling 
of the ghost-gluon vertex this fixed point can 
be explicitly calculated using propagator dressing 
functions alone. Employing a bare ghost-gluon 
vertex one obtains a^'^'iO) « 8.92/A^c M Re- 
cently it has been shown that this value together 
with the infrared exponent k « 0.595 are invari- 
ant in a class of transverse gauges that interpolate 
between Landau and Coulomb gauge 1131. 

2. Unquenching effects in propagators 

Dynamical chiral symmetry breaking is, be- 
sides confinement, the most important low energy 
property of QCD. It is a truly nonperturbative 
effect, since there is no dynamical mass genera- 
tion at every order in perturbation theory. On 
the quark level, the Schwinger-Dyson formalism 
is well suited to investigate the chiral symmetry 
breaking pattern also in the chiral limit of van- 
ishing current quark masses. It is therefore com- 
plementary to lattice simulations, which provide 
reliable results for large quark masses and vol- 
umes, but are yet severely restricted close to the 
chiral limit. In the SDE-framework the effects of 
dynamical chiral symmetry breaking in quenched 
and (partially) unquenched QCD have been in- 
vestigated in refs. |14ll5j . Based on the analytical 
results summarised in the last section, ansdtze for 
the vertices have been constructed such that the 




Figure 4. A diagrammatical representation of 
the coupled system of ghost, gluon and quark 
Schwinger-Dyson equations and the meson Bethe- 
Salpeter equation. Filled blobs denote dressed 
propagators and empty circles denote dressed ver- 
tex functions. 



system of SDEs for the ghost, gluon and quark 
propagators are closed and can be solved numer- 
ically (for details of the truncation scheme see 
refs. |14ll6p . Here we focus in particular on un- 
quenching effects in these propagators, which are 
generated by quarks in the gluon SDE, cf. fig. ^ 
Numerical solutions for the ghost and gluon 
propagators can be seen in fig. El In the infrared, 
the numerical SDE-results reproduce the analyti- 
cal power laws, eqs.Q. (This can be seen explic- 
itly on a log-log-plot, displayed e.g. in Ref. 
In the ultraviolet they reproduce the correct one- 
loop running from resummed perturbation the- 
ory. Compared to the results of recent lattice cal- 
culations 17 18 19^ (see also [2D]) we find good 
agreement for large and small momenta. Small 
deviations for the value of the infrared expo- 
nent K between continuum-SDE and lattice re- 
sults on a finite volunie may have methodical 
reasons. The resulting running coupling on the 
lattice does not reproduce the fixed point be- 
haviour in the continuum but vanishes in the in- 
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Figure 5. Comparison of the quenched and 
unquenched ghost and gluon dressing functions 
with recent lattice data |17I18I19] . The sea-quark 
masses are niu/d — 16 MeV, mg ~ 79 MeV in the 
lattice simulations and TOk/;; ~ 3.9 MeV, nis — 
84MeV in the SDE-approach. 



frared |1 712 1122] . These effects are also seen when 
one solves Schwinger-Dyson equations on a torus 
and are discussed in detail elsewhere [23 (see 
also [21). In the intermediate momentum region 
one clearly sees unquenching effects in the gluon 
dressing function due to the formation of quark- 
antiquark pairs from the vacuum. The screening 
effect from these pairs decreases the bump in the 
gluon dressing function considerably. The over- 
all difference in the size of the bump between the 



Figure 6. Comparison of the quenched and un- 
quenched quark mass (upper diagram) and wave- 
functions (lower diagram). 



lattice and the SDE-results is a measure of the 
influence of the (neglected) gluon-two-loop dia- 
grams in the gluon-SDE. 

The up/down and strange quark propagator 
functions, 

S{p)^Zf{p')/{-z^+M{p^)) (10) 

are plotted in fig. We clearly observe a large 
amount of dynamical mass generated in the in- 
frared. This mass is reduced by roughly 10 per- 
cent once quark- loops are taken into account. 
For large momenta the numerical solutions re- 
produce the logarithmic running known from re- 
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summed perturbation theory. There are notice- 
able unquenching effects in the intermediate mo- 
mentum region, which are, however, much smaller 
in size than those observed in the gluon propaga- 
tor. By explicitly solving the quark-SDE in the 
complex plane one finds a pair of complex conju- 
gate poles at (0.47 ± 0.29z)GeV (quenched) and 
(0.45±0.27i)GeV (unquenched) . Thus unquench- 
ing hardly has any effect on the position of these 
singularities (see however ref. 1^ for a discussion 
of the possible influence of scalar tensor pieces in 
the quark-gluon vertex on the analytical structure 
of the propagator) . 

3. Unquenching effects in light meson ob- 
servables 

The results for the gluon, ghost and quark 
propagators, discussed in the last section, serve as 
input for a calculation of light meson observables 
employing a Bethe-Salpeter equation (BSE), cf. 
fig. 0] The crucial link between the bound states 
and their quark and gluon constituents is pro- 
vided by the axialvector Ward-Takahashi iden- 
tity. It relates the quark self energy with the 
quark-quark interaction kernel in the BSE and 
thereby guarantees the Goldstone nature of the 
pions and kaons For details of the imple- 

mentation of this identity within the truncation 
scheme discussed here see ref. Results for 

the masses of the pseudoscalar and vector meson 
(in the isospin symmetric limit) as a function of 
the current quark mass are shown in fig. For 
ruq — > (the chiral limit) we observe a mass- 
less pion as expected. As for unquenching effects 
both the pseudoscalar and vector masses with 
larger quark masses are increased ^ 30MeV when 
quarks loops are taken into account. 

The values for the current quark masses at the 
physical point together with the resulting pseu- 
doscalar meson masses, leptonic decay constants 
and vector meson masses are given in Table ^ 
When fitted to the experimental pion and kaon 
masses the resulting up/down and strange-quark 
masses are lowered when quark loop effects are 
taken into account. This has also been observed 
in corresponding lattice simulations 27 28 . We 
furthermore see that the results for Jk and nip 
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Figure 7. Pseudoscalar (upper diagram) and vec- 
tor (lower diagram) meson masses as functions of 
the quark mass parameter. We compare results 
for the quenched and unquenched theory using 
three degenerate sea quarks. 



are quite insensitive to whether or not the sys- 
tem is unquenched (in the restricted sense of fig. 
4) . This leads to the conclusion that once the in- 
teraction has been fitted to the pseudoscalar ob- 
servables, the vector meson mass is largely fixed. 
The likely explanation for this is that the ground 
state pseudoscalar and vector mesons are both 
states with the lowest orbital angular momentum 
{L = 0) in the sense of the naive (quantum me- 
chanical) quark model ~ meaning that they are 
determined largely by the lowest spin contribu- 
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rus 


88.2 




86.0 


80-130 
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140.0 


139.6 
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130.9 


131.1 


131.0 


130.7 




266 


271 


271 




rriK 


494.5 




493.3 


493.7 


fK 


165.6 




169.5 


160.0 


nip 


708.0 


690.0 


695.2 


770.0 



Table 1 

Parameter sets and results for TOtt, /tt, "rriK, /k 
and nip for the quenched case {Nf = 0), the un- 
quenched case with three degenerate 'sea'-quarks 
[Nf = 3) and the physical quark configuration 
case {Nf = 2 -f 1) with two up/down quarks and 
one strange quark. The quark masses and the 
condensate have been determined using a large 
renormalisation point and subsequently evolved 
down to the scale /i = 2 GeV according to their 
one-loop running. All units are MeV. 



tions of the kernel in the Bethe-Salpeter equa- 
tion given by the ladder approximation (which is 
used here). The interaction then plays the same 
role in both channels, hence the similarity in re- 
sults. Note that the p meson calculated within 
the framework of the truncated Bethe-Salpeter 
equation here refers to a pure quark-antiquark 
meson with no allowed decay channel. A first step 
towards including the non-trivial decay width of 
the physical p meson in this formalism has been 
made in ref. [3^ (see also ref j31, for a calcula- 
tion of (jjrpp using quenched Bethe-Salpeter am- 
plitudes). Here the decay of p meson will lead 
to an additional shift of the p meson mass peak 
as evident from dispersion relations. As a con- 
sequence our present results should not directly 
compared to experiment. Although the p mass is 
slightly low the dispersive corrections via its two 
pion decay (and the inclusion of quark-gluon ver- 
tex corrections) might yield a satisfactory answer. 

A suitable quantity to compare results from the 
Green's functions framework with lattice data is 
the vector meson mass as a function of the pseu- 
doscalar meson mass. Since both quantities are 




Figure 8. Vector meson masses as a function 
of pseudoscalar meson masses. We compare the 
quenched and degenerate unquenched cases. The 
lattice results are taken from refs. |27l28j . 



physical no scheme ambiguities arise. The results 
are shown in fig.|Hl From the SDE/BSE-approach 
we compare results from the full SDE-setup, de- 
scribed above, with those obtained employing a 
phenomenological model for the quark-gluon in- 
teraction (the details of the model are described 
in ref. ^]). The model interaction and the effec- 
tive interaction of the full SDE-setup are comple- 
mentary to each other in the sense that the model 
interaction is confined to a quite narrow momen- 
tum region, whereas the interaction of the full 
SDE-setup has considerable strength in the in- 
frared and extends into the ultraviolet according 
to the correct one-loop scaling known from per- 
turbation theory. Together, both setups represent 
a measure for the theoretical error of our calcu- 
lation. This error is obviously of the same size as 
the combined systematic error of the different lat- 
tice simulations. In general, the results are in nice 
agreement with the lattice data. For pion masses 
below 240 MeV, where no lattice data are avail- 
able, the results show a nonlinear dependence of 
the vector meson mass on the pseudoscalar one. 
The effect of unquenching - when viewed as a 
function of the pseudoscalar meson mass - be- 
comes the same for both schemes: the vector me- 
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son mass is slightly increased when quark loops 
are taken into account. This trend is also seen 
in the lattice simulations |27l28j . where the ef- 
fect is even more pronounced. However, these un- 
quenching effects are small compared to the dif- 
ferences between both, the truncation schemes we 
employed and the systematic errors of the lattice 
results. 
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